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Abstract 

A first attempt to understand hadron dynamics at low energies in 
terms of the fundamental quark and gluon degrees of freedom incor- 
porates the effects of the gluonic field into a potential depending only 
on the spatial positions of the quarks, which are considered in the in- 
finite mass limit. A suitable framework for calculating such potentials 
between static quarks, i.e. a generalization of the Wilson loop will be 
discussed. 

Making a connection with recent Monte Carlo lattice simulations 
for the lowest two energies of a system of two quarks and two an- 
tiquarks, the static qqqq-potential will be calculated in perturbation 
theory to fourth order. The result will be shown to be exactly equal 
to the prediction of a straightforward two-body approach, which in 
Monte Carlo lattice simulations has been found to be a reasonable 
approximation for very small interquark distances. 
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1 Introduction 



As QCD is the accepted theory of the strong interactions, it is no doubt de- 
sirable to understand all hadronic phenomena directly in terms of the funda- 
mental fields of QCD. However, QCD being asymptotically free, perturbation 
theory is applicable only for very short distances and cannot cover the com- 
plete range of interest. At present lattice gauge simulations are the only way 
to study such systems. In a first approach, the static approximation is the 
natural choice, where the gluonic degrees of freedom are integrated out, and 
quark loops are ignored (the quenched approximation), giving rise to a po- 
tential between the stationary quarks. The potential of the quark- ant iquark 
system, where this approach — leading to the familiar Wilson Loop — is very 
well known, has been calculated extensively in Monte Carlo lattice simula- 
tions. (For recent data, see e.g. The ground state potential of this static 
system has also been calculated in perturbation theory upto sixth order ||. 

Here we shall describe how to generalize this procedure to multi-quark 
systems, especially to (gg) fc systems. However, even with present-day com- 
puters, qq lattice simulations are still very demanding, and the amount of 
computations needed increases rapidly with the number of interacting quarks. 
Reliable models for multi-quark systems expressing their potentials e.g. in 
terms of the well known gg-systems would therefore be of great help. Such 
two-body approximations have proven successful in many areas of physics, 
and these models can be formulated without difficulty. For the gggg-system, 
which is the simplest one that can be considered consisting of two colour sin- 
glets, this model has been tested against numerical data from a Monte Carlo 
simulation Q. For small distances the agreement has been found reasonable. 
It has also been observed ||] that the two-body model corresponds to lowest 
order perturbation theory. We shall be able to show that it is correct even 
to fourth order. To sixth order, however, three- and four-body forces begin 
to appear. 

2 The Generalized Wilson Loop 

While the concepts discussed below are of course well known in the context of 
the Wilson Loop for the gg-system, we find it useful to start with rephrasing 
these concepts in the case of an arbitrary number of quarks, leading to a 
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study of more complicated systems. 

When we have assembled a system of several quarks (and antiquarks), 
gluons will mediate a force between them. Treating this system in an ap- 
proximation as a quantum mechanical system of several static quarks, the 
interactions between the quarks are incorporated into a potential. This as- 
sembly of quarks is then expected to propagate in time with the usual factor 
of e~ ltH , where the interesting piece of the Hamilton operator H is the poten- 
tial energy. Thus, by calculating appropriate Green functions, the potentials 
of eigenstates of H can be extracted. 

2.1 Setting up Gauge Invariant States 

Because of confinement, it makes sense only to talk about systems of quarks 
where the overall states have colour singlet quantum numbers. The problem 
with setting up say a gg-system in a singlet is that the quark and antiquark 
are located a distance apart. This problem can be overcome by inserting the 
path ordered exponential U(x, y, A) = Ve 9 ^ T Aa ^ dZ)1 between the locations 
x and y of the quarks in the presence of the gauge potential A. Here g 
denotes the coupling constant and T a the representation matrices. Thus 
^(x)U(x,y)i(}(y)\0) will serve as a basis state in this case. We must also 
know how many basis states there are. When dealing with Green functions 
coming from Monte Carlo lattice simulations, they will have contributions 
from excited states of the gluonic field, and there are infinitely many of them 
even in the qq-ca.se. With suitable methods, the lowest potentials can be 
extracted, and several have been calculated for the quenched gg-system - see 
for example 0,0. The situation is different for Green functions calculated 
in perturbation theory. Here, unlike the lattice simulations, we can and must 
go to the infinite time limit. We do not expect to reach excited states of the 
gluon field in finite order perturbation theory, and thus the number of basis 
states for a system of several quarks is given by the usual arguments of group 
representation theory, e.g. one for the gg-system and two for the gggg-system. 
In the large time limit we expect that the effects of 'introducing' the quarks 
into the vacuum will be irrelevant in comparison to their time evolution, and 
the notion of a potential makes sense. 
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2.2 Diagonalization 



It may be shown that the state 



quarks and antiquarks qj at time t) 



4>i(-t,Xi) . . .U(-t,Xi,yi,A) . ..ip(-t,yi)\0) 



(1) 



satisfies Schrodinger's equation. Forming the overlap of states at time —t and 
t, we get an equation between Green functions and expressions of the form 
Aij d = (Ai\e~ ltH \Aj) , where \Ai) stands for some basis state and we have 
introduced the matrix A. By assuming a decomposition of these basis states 
into eigenstates of H, a diagonalization procedure will yield the potentials. In 
the case of the Green functions coming from lattice simulations, one considers 
a practical number of basis states, expands them in energy eigenstates and 
drops contributions with e~ ltEi for energies E^ above a certain limit. Of course 
we implicitly assume Wick-rotation. In perturbation theory, where a power 
expansion of e~' ltEi ^ in the coupling g will not be exponentially damped, 
we need to consider all linearly independent basis states, a number that is 
finite, as remarked in the last section. Because of this finiteness, we can find 
an invertible transformation to energy eigenstates, and the diagonalization 
is straightforward. In fact, given a matrix A satisfying certain consistency 
relations, we can perturbatively prove the existence of a time-independent 
basis transformation such that in this new basis A is not only diagonal, but 
its eigenvalues are of the form e~ ltEi<y9 \ Here the energy E^g) of the z-th 
basis state, which can be calculated perturbatively, is for static quarks equal 
to the i-th potential (apart from an irrelevant constant, the rest mass). 

2.3 Loops 

What remains to be done is to bring the Green functions of the last paragraph 
to more familiar forms. Since we work within the static approximation, the 
full quark propagator in the presence of gauge fields can be calculated ||: 




~ y ° ^5(x — y) x 
Q(y° - x )] 



i + i 
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Q(x°-y°) + 



(2) 
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We shall now outline how various contour integrations, i.e. loops arise. 
Considering the well-known gg-case, we find a path-ordered line integral from 



and the antiquark propagating backwards in time close the rectangle of the 
familiar Wilson loop. Starting with the Green functions described below 
eq. ([l|) and evaluating them for propagation from —t/2 to t/2, the following 
diagrammatic rule for calculating the Green function dealing with an arbi- 
trary number k/2 of quark-antiquark pairs (i.e. k quarks and antiquarks) 
partitioned into qq singlets is seen to hold: 

1. Draw two horizontal lines, the lower denoting time —t/2, the upper 
t/2. Mark the position of every quark and antiquark on the lower line 
and once again vertically above it on the upper line. 

2. At the —t/2 level connect every quark-antiquark pair that is set up as 
a singlet at —t/2 with a line, having an arrow pointing from antiquark 



3. At the t/2 level connect every quark-antiquark pair that is set up as 
a singlet at t/2 with a line, the arrow in which points from quark to 
antiquark. 

4. Join the quarks at the —t/2 level with quarks at the same position at 
the t/2 level, arrow pointing upwards, i.e. forward in time. 

5. Join the antiquarks at the t/2 level with the antiquarks at the —t/2 
level, arrow pointing downwards, i.e. backwards in time. 

6. Associate a path-ordered exponential of e %9 §c T A Z( z ) dz » together with a 
trace for every closed loop C occurring. 

7. Determine the overall sign: If the pairings at the —t/2 level are the 
same as those on the t/2 level, there must be a + sign. (This follows 
from the positivity of the norm on a Hilbert space if one lets t — > 0.) 
If this is not so, determine the sign of the permutation of antiquarks 
on the upper line that is necessary to give the same pairings as on the 
lower line. This is the overall sign. 




to quark. 
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8. Multiply by (5(0)e im *) fc , where k is the total number of quarks and 
antiquarks. 

9. Insert the factor so obtained in the numeratoifl of ^ DA ^ Dl]a ^ Dva ^ e — 

U f[DA % ][Dr l * a ][Dr la )e i J d4 * lC] 

This gives the Green function in the chosen singlet structure. 



3 The qqqq- Potentials 

In S77(N) gauge theory with quarks in the fundamental representation, we 
want to calculate the gggg-potential in perturbation theory to fourth order. 
It has been remarked in subsection [2.1| that there are two independent basis 
states for this system, and one easily recognises a choice of these in the 
two possible ways of pairing the system into two quark- ant iquark singlets. 
Assuming the first static quark at position Ri, the second at R2, and the 
antiquarks at R3 and -R4, we will label the two states \Ai) = I13I24 and 

1^2) = 3-141-23- 



3.1 Calculating the Green Functions 



According to subsection 2.3, we encounter the following types of loops: 



R \ 



R, 



R* 



t/2 



-t/2 



Ri 



C 



(Ai,-t/2|i4i,t/2) 



and 



1 With rj we denote the ghost fields, with C the Lagrangian without fermions 
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t/2 



(A 2 -t/2\A 1: t/2) 



-t/2 



Ra 



R* 



Ri 



In calculating the Green functions, we will adopt dimensional regular- 
ization in dimension D = 4 — 2e with a mass scale M. A wave-function 
renormalization will remove infinities associated with a diagram of the form 



- 3 



while a coupling constant renormalization is needed to make the sum 




T2 
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finite. The full gluon propagator is renormalized in the usual way. Note 
that, contrary to a smallest distance regularization frequently encountered 
in this context, there is no renormalization of the quark mass m. Making 

2 

the expansion in a — rather than in g and also expanding in time, a 
calculation yields for 

A = A(l) + aA(a) + a 2 A(a 2 ) + atA(at) + a 2 tA(a 2 t) + M!a(^^) (3) 



the expressions: 
A(l) = 



N 2 -N 
-N N 2 



A(a) 
A(aT) = 



N 2 {Un{M 2 R 13 R 2A )} —N {in ( ^^glj 4 ^ ) -ei-e 2 ] 

-ei-ea} TV 2 {41n(M 2 J R 14 i?2 3 )} 



AT 

_ f N 2 {-V sl } -N{-V sl - V s2 + V d } 

-N {-V Bl - V s2 + V d } N 2 {-V s2 } 



A(a 2 T) = 

/N 2 {-W sl ln(M 2 R 13 R 24 ) 
+ l~V s2 + V d ] x 

M«)-*]} 



-N{[-V sl -V s2 + V d ] x 

fin f^L^Ml) _ ei _ 62 " 

+ it s ^ Vd - y si)+ c 2(V d - Via)] 
- ^ [ln(i? 12 i? 34 )(V 5 i + V 52 -2y d )] 
+ ^[ln(i? 13 i? 24 )(V 52 -V i )] 
\ +^ 7 [ln(i? 14 i? 23 )(V 5l -V i )]} 



-^V{[-Vsl-Vs2 + V/]X \ 

"in f^Lpk^LMl) _ ei - e 2 l 

V K 12 K 34 / J 

+ f 7 [ei(y d -Vsi) + e 2 (V/-Vs2)] 
- fa [ln(i? 12j R 3 4)(Vsi+Vs2-2Vi)] 
+ ^ [ln(i? 13 i? 24 )(Vs2 - V d )] 

+ ^-[ln(i? 14 i? 23 )(Ki-Vd)]} 
^{-4^2 In (M 2 R 1A R 2Z ) 



/ 
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A (\{aTf 



( n2 {vA + ^(Vs2 ~V d ) 2 } -N{(V sl + V a2 -V d ) 2 + ^x 



\ 

[v d v sl + v d v s2 - v sl v s2 - vj] } 

-N{(V sl + V s2 - V d f + N 2 [V 2 2 + fa - V d f) 

V [VdVei + v d v s2 - v sl v s2 - vl] } J 
where, to save writing, the following shorthand notation has been adopted: 



edge(R 3 , Rx, R 4 , R 2 ) 



def f f 1 / 2 (R 3 -Ri)(R 4 -R 2 ) dwdx 



-1/2 



^^+w{R 3 -R 1 )-^--x(R 4 -R 2 ) 



2 



ei d = edge(R ± , R 3 , R 2 , R 4 ) 

e 2 = f edge{R u R 4 ,R 2 ,R 3 ) 

Vsi = V(R 13 )+V(R 24 ) 

V s2 = f V(R 14 )+V(R 23 ) 

V d = f V(R 12 )+V(R U ) 



ClS ao — J^cdfacdfbcd ~ N5 a b 

c 2 5 ab = Tr[T a T b ] = \5 ab 
c 3 l = E a T a T a = %ll 



2N 

and V(Rpq) is the two-body potential between a quark 'p' and an antiquark 
'q' a distance R pq apart. Diagonalization yields for the two possible energy 
eigenstates the two potentials correct to fourth order: 



(iV 2 -2) (V sl + V s2 ) + 2V d - jVy/iV 2 (V al -V a2 y + 4 (V al -V d )(Vs2-Vd 
" 2(iV 2 -l) 



(iV 2 -2) (K1 + K2) + 2Vd + iV y/iV 2 {V sl -Vs2) + 4 (V al -V d )(V S 2-Vd) 
1 2(iV 2 -l) 

(4) 

These potentials are exactly equal to those given by a naive two-body model 
(see 0) — and are one of the main conclusions of this work. 
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3.2 Consistency Relations and Edge Effects 



When integrating out the gluonic degrees of freedom and reducing a field the- 
ory to a quantum mechanical system, we have made assumptions which may 
need some justification. This is especially so in the perturbative case where 
we have postulated a finite number of energy eigenstates. We, therefore, 
want to have a look at the consistency of this approach. 

There are certain relations which must hold in order to guarantee proper 
exponentiation. We have remarked at the end of subsection |2]2| that after 
diagonalization we expect the diagonal entries of the matrix A to be (apart 
from a normalization factor) of the form e~ ltE ^ 9 \ An actual calculation will 
give a power expansion in g and t for these diagonal entries. This will not 
only determine the energy Ei(g) as the coefficient of t, but will also allow 
us to check the consistency of our calculation by inspecting the relations 
between the coefficients of higher powers of t. We were able to verify these 
consistency relations in our calculation to fourth order. 

Another important point to notice is that the terms we have abbreviated 
edge (R p , R q , R r , R s ) have cancelled in eq. (ffl). These terms are coming from 
the space-like contour integrations. In subsection |2. 1| we had introduced these 
space-like contour integrations via the path-ordered exponential U in order 
to establish colour singlets. The exact integration contour in U was of course 
arbitrary, and the physical potentials must not depend on it. Another way 
of looking at it is that in the large time limit the processes associated with 
bringing the quarks into their position in the distant past (and removing them 
again in the distant future) should become irrelevant. One may also make 
an adiabatic argument, turning on the coupling g in the distant past and off 
again in the distant future, thus making the notion of a colour singlet state 
at these times well-defined in spite of the spatial separation of the quarks. In 
any case, for the potentials of eq. (Q) to be meaningful, they must not involve 
terms originating in contour integrations in the distant past or future, and 
indeed they do not. 



3.3 Three- and Four-Body Forces in Sixth Order 

The fact that a straightforward two-body model is correct also to next-to- 
leading order may be surprising in light of the non-abelian nature of QCD. 
Hence we want to mention that we believe that the two-body model fails at 
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sixth order as three- and four-body forces appear at this order. While many 
diagrams that seem to give rise to such deviations from the two-body model 
actually vanish, we see no reason why for example effects from the following 
two diagrams should cancel in the calculation of the four-quark potentials to 
sixth order. 





where 



o 



stands for 




Hence three- and four-body forces will be introduced. In general, their nature 
seems to be complicated, but for some geometries simplifications are possible; 
e.g. for the four quarks on the corners of a regular tetrahedron there will 
be no contribution from quark self- interactions to four-body forces to sixth 
order. 
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4 Relation of our Result to Lattice Simula- 
tions 

Looking at the Monte Carlo lattice calculations for the gggg-system discussed 
in H , we observe that for small interquark distances of a few lattice spacings 
(with a lattice spacing of a ~ 0.12 fm) the two-body model gives a reason- 
able approximation in the sense that the four-quark potentials calculated 
from eq. (H) using the Monte Carlo two-body potentials are comparable to 
the four-quark potentials from the lattice simulation. The agreement im- 
proves the smaller the distances get. By comparing the perturbative (i.e. 
1/i?) and non-perturbative (i.e. linear) part in the usual parametrization of 
the gg-potential, one would expect to start entering the perturbative regime 
at distances of about two lattice spacings. However, at that stage the ap- 
proximation provided by the two-body model is already very good. The fact 
that the two-body model is correct to fourth order in perturbation theory 
certainly suggests that it should be a reasonable approximation in the per- 
turbative domain. So our result supports the belief that the results of the 
lattice simulations for small enough distances indeed are correlated to con- 
tinuum perturbation theory, and thus that continuum physics is extracted 
from the Monte Carlo calculations. 
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